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Abstract: We demonstrate that some of the suggested five supergravity points for study 
at the LHC could be approximately derived from perturbative string theories or M-theory, 
but that charge and colour breaking minima would result. As a pilot study, we then 
analyse a perturbative string model with non-universal soft masses that are optimised in 
order to avoid global charge and colour breaking minima. By combining measurements 
of up to six kinematic edges from squark decay chains with data from a new kinematic 
variable, designed to improve slepton mass measurements, we demonstrate that a typical 
LHC experiment will be able to determine squark, slepton and neutralino masses with 
an accuracy sufficient to permit an optimised model to be distinguished from a similar 
standard SUGRA point. The technique thus generalizes SUSY searches at the LHC. 
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1. Introduction 

The purpose of this work is to extend the discussion of LHC supersymmetry (SUSY) 
searches to include string models. We begin by discussing whether string models can be 
used to motivate previous work on LHC SUSY searches, and then suggest a well- motivated 
non-universal string model for a new pilot study. We go on to examine how the SUSY 
particles can be detected and how the model can be distinguished from a similar well 
studied supergravity (SUGRA) model. We reconstruct sparticle masses by looking for 
kinematic edges in qi X29 ~^ ^r^^Q ~^ Xi^^^^Q and ql — > X29 ~^ Xi^Q ~^ Xi^^^^Q decay 
chains, and in doing so generalize the method of |j^] by unifying the treatment of light and 
heavy sleptons. Additionally, with a novel method based on ||^], we further constrain the 
Xi and Ir masses by studying the kinematics of events containing pair produced sleptons: 
(gg/qq) — l^l^ l''^Xi^~Xi- In particular, this allows the mass difference between the 
Xi and Ir to be determined with sufficient accuracy to permit discrimination between the 
string model and the most similar standard SUGRA model. We suggest that our analysis 
is likely to be applicable, not just to string motivated non-universal models, but to other 
non-universal models as well. 
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2. Theory 



2.1 Introduction 

Throughout this work we assume that the effective theory describing TeV scale physics 
is the R-parity conserving minimal supersymmetric standard model (MSSM). Within this 
framework, the collider phenomenology is strongly affected by the SUSY breaking terms 

1 ^ 

/: = - ^ MaXala " J2 l'^*!^ " i^ijkWijk + Bl^HiH2 + H.C.), (2.1) 



a=l 



where i,j,k = Ql, d^, L/,, Hi, H2 and (pa, Aa are the scalar and gaugino fields of the 
MSSM (see e.g. Ref. |^). Wijk are the trilinear pieces of the MSSM superpotential, which 
written in terms of superfields is 

W = YeLHiE + YdQHId + YUQH2U + fiHiH2, (2.2) 

where we have suppressed all gauge and family indices. Y^, Y^;, Yu denote the down 
quark, charged lepton and up quark Yukawa matrices respectively. 



In Eq. (2.1), a general parameterisation of possible SUSY breaking effects has been 



employed^. Usually, the parameters are constrained by the condition of universality, 

rrii = mo, Aijj, = Aq, Ma = M1/2 (2.3) 



deriving from simple SUGRA models. Eq. ( |2.3| ) is subject to radiative corrections and 
should be imposed at the string scale M5. In the usual formulation of perturbative string 



theory, this corresponds to Ms ~ 5 x 10^"^ GeV, but Eq. (2^) is usually applied at the 
grand unified scale Mqut ~ 2 x 10^^ GeV as an approximation. Once B and /i are 
constrained by radiative electroweak symmetry breaking |^, the SUSY breaking sector is 
then characterised by one sign: sgn/i, and four scalar parameters: mo, -^1/2) and tan/3, 
the ratio V2/V1 of the two MSSM Higgs vacuum expectation values (VEVs). Once these 
are specified and current data are used to predict supersymmetric couplings such as the 
top Yukawa coupling and gauge couplings, the sparticle spectrum and decay chains are 
specified. Five points (denoted S1-S5) in the SUGRA parameter space mo, Aq, Mi/2) sgn//, 
tan (5 have been suggested for study of SUSY production at the LHC |^ and are catalogued 
in Table ||. These models have been weh studied in the context of the LHC @, |, ||, |l|], 
and we shall use them as a reference to compare and contrast with new models, which do 
not necessarily obey Eq. (|2.3|). 



2.2 SUGRA point compatibility with string models 



In Refs. |11, 12 1, the authors study the phenomenological viability of string and M-theory 
scenarios coming from the desirable absence of dangerous charge and colour breaking (CCB) 
minima or unbounded from below (UFB) directions in the effective potential. One of the 

^It has been assumed that non-standard terms such as those discussed in Ref. are disallowed because 
they cause a naturalness problem in the presence of gauge singlets. 
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Table 1: Compatibility of LHC SUGRA points S1-S5 with strongly /weakly coupled string models. 
A tick in the 'Weak'/'M-theory' column indicates that the sparticle spectrum could approximately 
be derived from the weakly/strongly coupled string models respectively. 



models considered in [12| is weakly coupled string theory with orbifold compactifications. 
In this case, the soft masses evaluated at the string scale are dependent upon the modular 
weights rii of the fields and do not necessarily conform with Eq. (2.3) Other string 
scenarios exist [14, 15, 16 1 which erase the UFB/CCB global minima which we do not 
explicitly investigate. For all modular weights equal to -1 however, the tree-level pattern 
of soft-masses conforms with Eq. (2.3), with the additional constraint 



M- 



1/2 



-Aq = VSruQ. 



(2.4) 



In Table |^, we display under the "Weak" column whether each standard SUGRA point is 
approximately compatible with this sub-class of universal string models. None of S1-S5 fit 
exactly with this scenario, but S5 is the closest and would have a similar sparticle spectrum 
if TTT-o were 173 GeV instead of 100 GeV. In regions of consistent radiative electroweak sym- 



metry breaking (REWSB), the string model version of S5 has dangerous UFB minima |12|. 
We reject this class of model for further study, partly because S5 gives a similar spectrum, 
and partly because it is already well studied, but mainly because of the UFB problem in 
the potential mentioned above. 

The "M-theory" column of Table |l| refers to compatibility with the strong-coupling 
limit of Eg X Eg Heterotic string theory |17]. Given simple assumptions (that SUSY is 



spontaneously broken by the auxiliary components of the bulk moduli superfields), the soft 
terms of M-theory valid at Mqut are [18| 
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(e(6 + e) sin^ 6* + (3 + 2e) cos^ 6 - 2V3e sin 6 cos I 



3/2 



3 + e 



(3 - 2e) sin 9 + VSe cos a 



(2.5) 
(2.6) 

(2.7) 



So, the SUGRA parameters M1/2, tuq, Aq become replaced by the goldstino angle 
sin^, the gravitino mass M3/2 and a ratio of moduli VEVs 0<e<l. e = corresponds 
to the weakly-coupled perturbative string. For S1-S4, we notice from Table H that Aq = 0, 
allowing us to solve Eq. (2.7): 
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tan 6 



(21 



2e-3 

Specifying e then determines the ratio |mo/Afx/2p) which is displayed in Fig. |^a. 
|mo/Mi/2p is always greater than one except as e ^ cxd, ruling out M-theory derivations 
of SI and S2. S3 and S4 are compatible with e = 0.21 and 0.53 respectively, as indicated 
in the figure. 
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Figure 1: Spectrum of M-theory valid at Mqut- (a) is valid for S1-S4, whereas (b) is valid for 
S5. Dotted lines show the relevant parameters for models that reproduce the LHC SUGRA point 
spectrum (S3-S5). 



Eq. (2.8) is not relevant for S5 because 7^ 0, but the condition M1/2 = ^0 may be 



solved yielding 



-e(3 + 2e) 



tan6' = ^ (2.9) 

V3(3 + e-e2) ^ ' 

I mo /Ml/2 P is then plotted using this relation against e in Figure lb. The fi gure indicates 
that M-theory does not reproduce S5 for any realistic value of e. 

To summarise, Table |l| shows that S3-S4 are compatible with Eq.s (2^)- (2^), and 
therefore M-theory. However, Refs. [0, 12 1 show that each of the points corresponding 



to S3-S4 is in conflict either with UFB, CCB or REWSB constraints. In fact, all of the 



M-theory parameter space examined in Ref. |12| was shown to be in conflict with one of 
these constraints. 

We have therefore shown that while the LHC SUGRA points include models which may 
be derived from weakly or strongly coupled strings, they possess potentially catastrophic 
global CCB and/or UFB minima. The existence of a global CCB or UFB minimum does 
not necessarily rule out a model. Some models have meta-stable minima with lifetimes 



longer than the current age of the universe |19, 2C, 21 1. The question of which minimum 
the VEV of scalar fields rest in is one of cosmology, and beyond the scope of this paper. 
We therefore take the view that if the bounds from CCB/UFB global minima are not valid, 
examples of weakly/strongly coupled string models are approximated by or included within 
S1-S5 so there is no need to construct another similar model to examine the string-derived 
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SUSY phenomenology at the LHC. If one should take the global CCB/UFB minimum 



bounds strictly however, Ref. |12] indicates that no variation of parameters in the class of 
string models considered above will result in a model without CCB/UFB problems. We 
thus conclude that it is useful to examine a new class of string model which does not have 
difficulty evading CCB/UFB constraints. We note that it is also possible to evade the 
CCB/UFB constraints by lowering the string scale in type I string models |15, 16|. 



2.3 Optimized string model 

We now turn to the analysis of a model specifically designed to provide a large region of 
parameter space without CCB/UFB problems [l^. It is a weak coupling model, where the 
modular weights have been chosen so that 

"Qi = "-d^ = nu<^j^ = -2, = "-e^ = riH^,2 = -1- (2-10) 

The model is non-universal, but still in a family independent way, and so avoids seri- 
ous problems associated with flavour changing neutral currents. As a consequence of the 
assignments in Eq. ( p. 10 ), the string scale soft masses are 



"^Hi,2 = ™L = "^e^ = ^3/2 sin^ ^> 
^LiHieJj = -Ml/2 = \/3M3/2sin6l, 

'^Ql = T^d-^ = "T-u^ = 0, At = M3/2(\/2 - — ), (2.11) 

so that the squarks are light and the sleptons heavy at the unification scale. This has the 
effect of ameliorating the CCB/UFB problem. It should be noted |12] that the A-terms 



of squarks other than the stop are not calculable for small tan/?. We shall approximate 
them here to be equal to At, but in fact they have a negligible effect upon the phenomenol- 
ogy/spectrum unless tan/? is large, in which case Ah = Af. 

To be definite, we choose model parameters tan/3 = 10, M3/2 = 250 GeV and 9 = 7r/4. 
We call this optimized model 01 and numerically its parameters are 



mHi^2 = = "T-e^ = 177 GeV, 
-AL^H.e'^^ = Ml/2 = 306 GeV, 

mq^ = md-=^ = = 0, Aq^h^u-^^ = ^QLHid-^ = 137 GeV. (2.12) 
We make the approximation that these relations hold at Mqut 

~ 2 X lO^*^ GeV, but 

it should be borne in mind that logarithmic corrections from renormalisation between 
-^string ~ 5 X 10^^ and Mqut are expected. The spectrum and decay chains of the 
sparticles are calculated using ISASUGRA and ISASUSY |22| using Eq. (|2T^ ) as input. 



Using three different Monte Carlo programs and the CTEQ3L parton distribution func- 
tions [|^], we calculate the total cross sections of SUSY particles at the LHC for models 
S1-S5 and 01. Table shows the comparison of HERWIG6 . 2 total cross-section with those 



calculated by ISAJET7.40 ||2|] and SPYTHIA ||5|. In each case, we have used the CTEQ3L 
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Table 2: Comparison of LHC SUSY production hard subprocess total cross sections (in picobarns) 
for the LHC SUGRA points S1-S5 and the optimized string model 01. By "associated production" 
we mean the production of a chargino or neutralino in association with a gluino or squark. CTEQ3L 
parton distributions were used and the statistical fractional error on each result is 2%. 

parton distribution functions and the spectrum is calculated by ISASUGRA with mt = 175 
GeV. As can be seen from the table, the three Monte-Carlo programs agree to about 10%. 

Fig. ^ displays the result of the calculation using HERWIG6.0. It is noticeable from the 
figure that S5 and 01 have broadly similar SUSY production cross-sections except for the 
sleptons which are noticeably lower. This can be understood by comparing the spectra of 
the two models, displayed in Table ^. The spectra are approximately similar, except for 
the sleptons which are heavier in 01 relative to the squarks. This is a consequence of the 
different choice of modular weights for the sleptons and the squarks. We therefore propose 
to analyse LHC SUSY production in 01 using S5 as a benchmark or comparison. We will 
show that the two models can be distinguished experimentally. 



^This version of HERWIG, not officially released, was a developmental version of HERWIG6. 1 
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Figure 2: LHC SUSY production hard subprocess total cross sections for the LHC SUGRA points 
S1-S5 and the optimized string model 01. HERWIG6.0 was used with the CTEQ3L parton distribu- 
tions. Statistical fractional error on each point is 2%. Comparison with SPYTHIA, ISAJET7.40 is 
shown in Table |[ 
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Table 3: Comparison of sparticle and Higgs spectrum of 01 and S5 (in bold type). These spectra 
were calculated using ISASUGRA and all masses are quoted in GeV. Sign conventions are as in 
ISAJET. The masses of the second family of sparticles are approximately equal to those of the first. 
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3. Experimental observability 
3.1 Method 

The primary experimental aim is to take previously developed model-independent methods 
for measuring SUSY particle masses (which were developed at standard SUGRA points) 
and by testing them in the context of a new optimised model, identify where these methods 
need to be generalised to perform well in both models. Secondarily it is to be shown 
that, after modifying these methods' treatments of the slepton sector, their performance is 
sufficient to distinguish between the optimised and standard SUGRA scenarios. 




Figure 3: "Sequential" decay. 

To accomplish these aims, some model-dependent assumptions have to be made. In this 
analysis, R-parity is taken to be conserved and certain sparticle decay chains are assumed 
to exist. In R-parity conserving (RFC) models, SUSY particles are only produced in 
pairs, and the lightest SUSY particle (LSP) is stable. SUSY events contain two of these 
LSPs which, being only weakly interacting, escape detection and lead to SUSY events 
with large amounts of missing transverse energy - the standard signature for RFC models. 
The fact that these massive particles go missing from all RFC SUSY events means that 
in these models it is not usually possible to measure particle masses by reconstructing 
entire decay chains. Thus, as in endpoints in kinematic variables constructed from 
SUSY decay chains must instead be examined. Specifically, the "sequential" decay mode 
QL X2I ~^ ^j?^near9 ~^ Xi^iaJnearQ (see Figure |^ and the "branched" decay modes 
ql X2I ~^ XiA"? ~^ Xi^^f'^Q (see Figure form the starting point for this investigation. 




Figure 4: "Branched" decays through the Z and Higgs bosons. 

The kinematic edges used in [|l|] to identify particle masses at S5 contain: 

• l^l^ edge: This picks out, from the "sequential" decays, the position of the very 
sharp edge in the dilepton invariant mass spectrum caused by X2 ~^ followed by 

i^ixl- 
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• l^l^q edge: In "sequential" decays, the llq invariant mass spectrum contains a linearly 
vanishing upper edge due to successive two-body decays. Our theoretical model of 
the l^l~q edge is more model independent than that used in earlier studies as it does 
do not assume any relation between the sparticle masses, other than the hierarchy 
< m^o < m^^ < m^o < rriq necessary for the "seqential" decay chain to exist at all. 

• l^l^q threshold: This is the non-zero minimum in the "sequential" llq invariant mass 
spectrum, for the subset of events in which the angle between the two leptons (in 
the centre of mass frame of the slepton) is greater than 7r/2. This translates into the 
direct cut on mu described in Section ^.4.1 . 

• hq edge: This is one of the two instances of the Xq edge. In general, the Xq edge is 
the upper edge of the distribution of the invariant mass of three visible particles in 
the "branched" decays of Figure ^ The position of this edge is again determined by 
two-body kinematics. Depending on the Higgs mass and the mass difference between 
the X2 2ind the Xij one of the two "branched" decay chains will be strongly suppressed 
with respect to the other, so typically only one edge will be visible. At S5 the Higgs 
mode dominates, while the Z mode will dominate at 01. 

As detailed later in this analysis, new edges are added to the list above. Two of them 
are more general versions of other edges treated in Q , while another is entirely new. They 
are the l^q high-, l^q low- and Mt2 edges respectively. To summarise, the overall method 
we apply in this paper consists of: 

• finding a model-independent set of cuts which can be used to select events from which 
the endpoints of all the kinematic edges may be measured, 

• obtaining an estimate of the accuracy with which the edge positions might be deter- 
mined by an LHC experiment, 

• performing chi-squared fits of the expected positions of these edges (as functions of 
the sparticle masses) to a set of "simulated edge measurements" one might expect 
from an ensemble of such experiments, and 

• interpreting the results as the statistics contribution to model independent sparticle 
mass measurements of a typical LHC experiment. 

Were all squarks to have the same mass, then (with perfect detector resolution and 
with infinite statistics) the endpoints of the edges listed above would be found at the 
positions given in Table ^. Note that these positions depend on four unknown parameters, 
namely the masses of the squark, the slepton and the two neutralinos participating in each 
decay. The mass of the lightest Higgs boson also appears, but we assume this will already 
be known or will be obtained by other methods (e.g. ||2^) to within 2%. 

In a more realistic "non-degenerate squark masses" scenario, every distribution with a 
squark related edge will actually be a superposition of many underlying distributions (one 
for each squark mass), each having its edge at a slightly different location. This results 
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in some smearing of the edges, even before detector effects (resolutions, acceptances, jet 
energy calibrations etc.) are taken into account. 
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Table 4: The absolute kinematic endpoints of invariant mass quantities formed from decay chains 
of the types mentioned in the text for known particle masses. The following shorthand notation has 
been used: x = 't^|o , I = "m-j , £, = 't^|o , <j = rn?q and X is m^^ or depending on which particle 
participates in the "branched" decay. 



It is not possible to measure each of the squark masses separately. Consequently the 
quantity rriq in Table ^ will, after being obtained from the "smeared" edges, represent a 
squark mass scale rather than a specific squark mass. In all squark related edges, except 
the l~^l~q threshold, the contribution to the outermost part of each edge is provided by 
the heaviest squarks. In the case of the l~^l~q threshold, the "true" endpoint is set by 
the lightest squark. However, if (as at 01 and S5) the other eleven squarks are much 
heavier (see Table ^), it is easier to measure the contribution coming from them. A strong 
correlation between rriq and the mass of the heaviest squark would therefore be expected. 

More work would be required to fully understand the "theoretical" systematic errors 
including the use of a full simulation of the expected edge positions. Better edge models 
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than those used later in this analysis will definitely be needed to permit measurements of 
edge positions to be better correlated with functions of the particle masses. Other sources 
of systematic errors which require further analysis are the detector effects mentioned above, 
combinatorial backgrounds near the edges and possible cut biases. Such systematic errors 
can only meaningfully be studied when real data are available. 

Although they are beyond the scope of this paper, it is assumed that the above inves- 
tigations could be performed so as to leave the eventual edge resolutions determined only 
by statistics and detector resolution. In this analysis, then, all edges are fitted with simple 
shapes (see Section 3.5) with the intention of extracting only an estimate of the uncertainty 



on the edge position, and not to obtain the edge position itself. 

The first four edges listed in Table ^ thus constitute a minimal constraint on the four 
unknown sparticle mass parameters, which may then be further constrained by the other 
measurements which follow. 

3.2 Other measurements 

3.2.1 Near, far, high and low l^q edges 

In "sequential" decays there are three observable outgoing particles. There are only four 
different ways of grouping these particles together, and so only four different invariant 
masses may be formed from them: mu, mug, mi^^^^g and mi^^^q. The first two of these, mn 
and raiiq^ may be formed from the observed momenta without knowledge of which lepton 
was ^near and which was /far; only the total lepton four- momentum is required. Edges in 
these invariant mass distributions have already been described. If it were possible to tag the 
near- and far- leptons separately, the third and fourth invariant mass combinations could 
also be formed unambiguously, allowing the positions of two more edges, w-™^^^ and ?Ti™'^^, 
to play a part in the final fit. However, such tagging is impossible. If further information 
is to be gathered from these decays in a model independent way, it is then necessary to 
look for edges in variables (functions of nii^^^^^g and nni^^^q) which are observable. On an 
event-by-event basis, then, we define 

"izq(high) = max {rrii+q, rui^q) = max {mi^^^^q, mi^^^q) (3.1) 

and 

"i/9{low) = min {mi+q, nii-q) = min (m/,^^^^g, mi^^^q). (3.2) 

The simplest theoretical predictions for the positions of the corresponding edges, T^Yq(high) 
and w-fqQo^), are listed in Table ^, along with the positions of the near- and far-edges. 
Note that although the position of the high-edge is just the higher of and rnf^^^, 

the position of the low-edge has a more interesting form. This asymmetry arises because 
it is not always kinematically possible for the invariant mass of the lepton/quark pair 
coming from the lower near/far distribution to approach min f^^^^qimf^^^ arbitrarily 
closely, while simultaneously requiring that this invariant mass is less than that of the 
other lepton/quark pair. 
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3.2.2 Zq edge 

Since the neutralino mass difference at 01 (109 GeV) is too small to permit X2 ~^ ^Xi> 
"branched" decays are mediated by the Z. (Particle masses may be seen in Table ^.) The 
cuts developed in [|l| for picking out this special case of the Xq edge at S2 are found to 
perform equally well at 01, so they are adopted unchanged. Although it might be possible 
to benefit from adapting these cuts to 01 slightly, this temptation is resisted in order to 
retain model independence. 

3.2.3 Mt2 

In order to constrain the slepton and neutralino masses better, we construct another vari- 
able, Mt2{x)i whose distribution relates just these two masses. We base this on the 
variable, proposed in [^, which looks at events containing two identical two body decays: 
^ X — > Yaia2 YbiCib2C2, where the particles of types a and c are of unknown 
mass, where particles of type c are undetectable, where the longitudinal momentum of the 
incoming particles is also unknown, and where it is assumed that Y does not contain any 
unobservable particles such as neutrinos. In such cases, the variable provides a kinematic 
constraint on the masses of a and c. We seek to apply this primarily to LHC dislepton 
events of the form qq — > /j^/^ ^^Xi^^Xi and so define our Mt2{x) by: 



where 



^T2(x) = min [max{mT.(p^,pipx),"Zj.(p^,pi2,x)}J (3.3) 
Pi+P2=Pt 



mUPT, qr, X) = mf + + 2(^^E^ - • qr), (3.4) 



E^j, = -^/p^^ + mf and = ,J + x^- (3.5) 

This definition includes the lepton masses for completeness, although these are neglected 
in all computations. 

The value MT2ix) takes for a given candidate dislepton event is a function of: the 
transverse missing- momentum vector, pij,; the transverse momentum vectors of the two 
leptons, p^ and pf^; and one other parameter - an estimate of the neutralino mass, x (iiot 
to be confused with the actual mass of the neutralino, m-o). Unlike the other parameters, 
the value of x is not measured in each event ~ events may be reinterpreted for different 
values of x- ^T2{x) has the property that, for signal events in a perfect detector, 

max MT2{'niyo) = "i,-. (3.6) 

events L Xi J 

Thus when y is indeed m^o, then the distribution of 11^2 has an end point at the slepton 
mass. Since the other observables allow m^o to be measured in a model independent 
way, Mt2 can then be included in the analysis to provide an additional constraint on the 
slepton/neutralino mass difference. 
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In practice the Mt2 edge can be distorted by the finite resolution of the detector 
or missing energy from soft underlying events. Additionally, standard model (SM) back- 
grounds provide constraints on minimum detectable slepton/neutralino mass differences 
(see Section ^^). However the most important factor to control is the ability to correctly 
identify which particle species is contributing to an observed Mt2 edge ~ particularly since 
it is possible to have multiple edges in the non standard model contributions to Mt2 distri- 
butions. At S5, for example, where both the right- and left-sleptons are lighter than at 01 
(see Table both Z^Z]^ and /^/^ events pass the cuts, and two edges are generated. The 
edge coming from the (lighter) right-slepton falls well under the SM background and so is 
not measurable, while the (heavier) left-slepton still has a cross section for pair-production 
high enough to let it form a good edge of its own at MT2(m^o) = w-^^- It is important 
that this edge is not mistaken for the Ir, so methods are needed to dismiss it. A detailed 
prescription of how to go about such a dismissal is beyond the scope of this paper, but it 
would clearly be accomplished by looking for inconsistency between a given edge-particle 
hypothesis and all the other sparticle masses, the branching ratios and (in particular) the 
strongly mass dependent pair-production cross sections, which could all be measured by 
other means. 

3.3 Event generation and detector simulation 

All events, except those from qq — > W^W~ background processes, are simulated by 
HERWIG6 . 0. The VF-pair events are generated by ISA JET7 . 42 ||22| . The detector chosen for 
simulation is the ATLAS detector ||2^, ^ ^ , one of the two general purpose experiments 
scheduled for the LHC. The LHC is expected to start running at a luminosity of about 
10^^ cm~^s~^ and this is expected to be increased over a period of about three years to 
the design luminosity of 10^"^ cm~^s~^. These two periods are referred to, respectively, as 
the periods of low and high luminosity running. 

The performance of the detector is simulated by ATLFAST2 . 16 |^0[ which is primarily 
a fast calorimeter simulation which parametrises detector resolution and energy smearing 
and identifies jets and isolated leptons, in both the low and high luminosity environments. 
Throughout this analysis, the parameters controlling ATLFAST's jet and lepton isolation 
criteria are left with the default values appropriate to the apparatus: i.e. jets must satisfy 
> 15 GeV and must lie in the pseudo-rapidity range — 5 < r/-' < 5, while electrons must 
have pip > 5 GeV, muons > 6 GeV and both must lie in —2.5 < if < 2.5. For lepton 
isolation a maximum energy of 10 GeV may be deposited in a cone about the lepton of 
radius 0.2 in (r/, (;/))-space {(j) being the azimuthal angle) while its separation from other jets 
must be at least 0.4 in the same units. 

At high luminosity approximately 20 minimum bias events ("pile-up" events) are ex- 
pected to occur in each bunch crossing. Pile-up events are not simulated by ATLFAST2. 16, 
although it does alter its reconstruction resolutions to reflect the two different luminos- 
ity environments. It must be checked that any cuts applied at high luminosity will not 
be affected substantially by pile-up events. Within this article, events corresponding to 
100 fb~^ are generated and are reconstructed in the high luminosity environment. This 
approximately corresponds to one year of high luminosity running. 
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3.4 Event selection cuts 



Section 3.4.1 summarises all the cuts used to obtain the edges. The l^l^ edge, l^l^q edge 
and l~^l~q threshold cuts do not differ significantly from those used in ||l|. The l^q edge cuts, 
however, do. The most significant change is the relaxation of the splitting requirement. 
The original splitting requirement tries to guarantee that the jet which comes from the 
quark produced in association with the observed dilepton pair is correctly identified. It 
achieves this by insisting that: 

• both the dilepton pair and one of the two hardest jets ji and j2 (ranked by pr) are 
consistent with being the decay products of a squark (i.e. rriuj. < mcutofr)i and 

• the invariant mass rriuj. of the two leptons and the other of the two highest pT jets is 
inconsistent with these being the decay products of a squark (i.e. muj. > mcutoflf)- 

Although the above demand for inconsistency increases the purity of the signal, it has a 
significant detrimental effect on the efficiency. In our analysis the consistency requirement 
is retained but the inconsistency requirement dropped. Instead we require that mu is 
inside the expected region, given the edge measurement. We also perform background 
subtraction, modelling the opposite-sign same-lepton-family (OSSF) backgrounds by the 
distributions produced by those opposite-sign different-lepton-family (OSDF) events which 
pass the same cuts. 

The production cross section for dislepton events suitable for analysis is, in all 
models, by far the smallest (see Figure §). So to show that Mt2 events passing cuts are 
not in danger of being swamped by small variations in backgrounds or the cuts themselves, 
two set of cuts referred to as "hard" and "soft" are developed. The "hard" cuts attempt 
to maximise the signal to background ratio in the vicinity of the dislepton edge in the Mt2 
spectrum, while the much looser "soft" cuts try to maximise statistics at the edge at the 
expense of allowing in additional SUSY backgrounds. 



3.4.1 Cuts listed by observable 

This section summarises the cuts used in the analysis. For notational purposes, recon- 
structed leptons and jets are sorted by pT- For example j2, with corresponding transverse 
momentum p^?, is the reconstructed jet with the second highest px- All cuts are require- 
ments unless stated otherwise. 



l^l edge 

'^leptons = 2, both leptons OSSF and p^r^ > p'''} > 10 GeV. 
njets > 2 and p^ > p!^ > 150 GeV. > 300 GeV. 

The kinematics of OSSF leptons coming from background processes which produce 
uncorrelated leptons (for example tau decays) are modelled well by OSDF lepton 
combinations. Consequently, edge resolution is improved by "flavour subtracting" 
OSDF event distributions from OSSF event distributions. 
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q edge 

'^leptons = 2, both leptons OSSF and 

Pt>Pt> 10 GeV. 
^^jcts > 4, > 100 GcV, and p^^ > >p'^> p^ > 50 GeV. 

> max (100 GeV, 0.2Meifective) and MeSective > 400 GeV, where Meffective is defined 
by the scalar sum: 

Meflfective = + P^ + P^ + P^ + p't ■ (3.7) 

Since the desired edge is a maximum, only the smaller of the two muq combinations 
which can be formed using ji or j2 is used. 

Flavour subtraction is employed here as at the edge. 

l^l~q threshold 

All the cuts for the edge as above. 

In addition m™^/\/2 < mu < m'^^, where the value of m™**^ would be obtained 
from the Z+Z~ edge in practice, but was determined theoretically in this analysis. 

Since the desired edge is a minimum, only the larger of the two mug combinations 
which can be formed using ji or j2 is used. 

hq edge 

'^leptons = and > 300 GeV. 

Exactly two b jets with > > 50 GeV. No other b jets, regardless of pT- 

At least two non-6 jets with > > 100 GeV, with at least one inside —2.0 < 
r]^ < 2.0. 

within 17 GeV of Higgs peak in m^j spectrum. 

Since the desired edge is a maximum, the non-6 jet (i.e. or jg^) chosen to form 
the nihq invariant mass is that which minimises nihg. 

Zq edge 

"leptons = 2, both leptons OSSF and > Pt > 10 GeV. 

mil within 2.5 GeV of the centre of the Z-mass peak in the mu spectrum. 

At least two non-6 jets exist with pj, > 100 GeV and '^t > 300 GeV. 

Since the desired edge is a "maximum", the jet chosen to form the mzq invariant 
mass is the one (from those with pj, > 100 GeV) which minimises mzg- 

Again, flavour subtraction is used. 

I'^q high and I'^q low edges 

All the cuts for the l~^l^q edge above are required. 

Additionally, events consistent with the edge measurement are selected by asking 
for mil < + 1 GeV. 
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To choose the jet from which to form mi±g we select ji = ji or j2, whichever gives 
the smaller value of mnj, and require mnj. < rricutoSi where mcutoff is chosen to 
be above the l~^l~q edge, but is otherwise arbitrary. For easy comparison with 
"^cutoflf = 600 GeV was used in this analysis. 

Finally, the two invariant mass combinations, ^i^;±g, are assigned to ?TT.zq(high) s-nd 
"T'«g(iow) as defined earlier by Equations (|3.1| ) and ( |3.2| ). 

Mt2 edge (hard cuts) 

Events are required to have exactly one OSSF pair of isolated leptons satisfying 
> 50 GeV and > 30 GeV. 

6t < 20 GeV is required, where 6t = IPj- + Pt + J^rl- Large 6t in signal events is 
indicative of an unidentifiable transverse boost to the centre of mass frame, perhaps 
due to initial state radiation. Large 6t in other (not necessarily signal) events simply 
suggests an inconsistency with the desired event topology. 

Events containing one or more jets with p^j, > 40 to 50 GeV are vetoed. This cut 
also helps to reduce standard model backgrounds, notably ti. The lower this cut is 
placed, the better for the background rejection. However the cut cannot be placed 
too low (especially at high luminosity) due to the significant number of jets coming 
from the underlying event and other minimum bias events in the same bunch crossing. 
For order 25 minimum bias events per bunch crossing, |31| estimates that about 10% 



(1%) of bunch crossings will include a jet from the underlying event with a pT greater 
than 40 GeV (50 GeV). Our results are not sensitive to variation of the jet veto cut 
between 40 and 50 GeV, where at least 90% signal efficiency is expected. 

Events with \mij^i^ — mz\ < 5 GeV are vetoed to exclude lepton pairs from Z bosons. 

m^jij > 80 GeV and > 80 GeV are also required. 

Mt2 edge (soft cuts) These are as above, but 

• the "^T requirement is lowered from 80 to 50 GeV, 

• the upper limit for St is extended from 20 to 90 GeV, and 

• the dilepton invariant mass cut is removed altogether. 
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3.5 Edge resolutions 




400 600 800 1000 
(c1) Highm, (GeV) 





Si 400 



400 SOO 800 1000 
(b) m,„ (GeV) 




400 SOO 800 1000 
(c2) Low m, (GeV) 




400 SOO 800 1000 
(e) m,„ (GeV) 



Figure 5: The S5 distributions wliose endpoints are described in Table ^: a) tiie edge, b) tiie 
l^l^q edge, cl) tire l^q high-edge c2) the l^q low-edge, d) the l^l^q threshold and e) the hq edge. 



Plots were produced with the cuts described in Table 3.4.1. The number of events corresponds to 
100 fb^^ of high luminosity running. 



By applying at S5 the cuts listed in Section 3.4.1| we reproduce the nin and mug 
distributions whose edges are analyzed in (see Figure |5|: a, b and d). We also confirm 
that the same cuts may also be used to generate edges of similar quality at 01 (see Figure 
|6|: a, b and d). In addition, plots cl and c2 in Figures ^ and |6| display the m/g(high) 
and m.ig(iow) distributions generated from the modified l^q edge cuts also listed above. 
Numerical results are summarised in Table ^. 
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Figure 6: The 01 distributions wliosc cndpoints are described in Table ||: a) the edge, b) the 
l^l~q edge, cl) the l^q high-edge c2) the l^q low-edge, d) the l^l^q threshold and e) the Zq edge. 
Plots were produced with the cuts described in Table |3.4.1 . The number of events corresponds to 
100 fb~^ of high luminosity running. 



Similar pictures are seen at S5 and 01. The most obvious difference between the two 
sets of data is the peak at the Z mass present in Figure ^a but absent in Figure This 
peak comes from direct X2 Xi-^ ~^ decays, which have a branching fraction of 

39% at 01 compared with 0.65% at S5. The most common direct neutralino decay mode 
at S5 is through the /iq ( 65% ) of which only a tiny proportion ( 0.02% ) decay to the two 
light leptons ~ the partial decay width r(/io — > //) going approximately as m^. Scenarios 
in which the edge happens to coincide with the Z peak require special treatment and 
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(a) (b) (c) 

Figure 7: Line-shapes used for fits to the mass distributions. Shape (a) (or its reflection in the 
vertical axis) is the standard "straight line" used to extract the l^l~q edge, l^q high-edge, l^q 
low-edge and l^l~q threshold resolutions. Shape (b) models the expected shape of the sharp Z+Z^ 
edge in the absence of detector effects. The Z+Z" edge is actually fitted with shape (c), which is 
identical to (b) except that it is smeared with a gaussian resolution whose width is a free parameter 
of the fit. 



are not considered here. 

In order to make useful statements about the degree to which model parameters may be 
extracted from the endpoints and edges of these distributions, it is necessary to obtain an 
estimate of the accuracy with which these observables may be measured. It is expected that 
the errors on all of the observables considered here will eventually be statistics dominated, 
so simple fits have been made to the data to obtain estimates of the statistical errors on the 
edge or end point locations. The shapes fitted to the data (see Figure 0) and the algorithms 
for determining the boundaries of the fitted regions have been kept as simple and generic 
as possible, with the intention of making the fit results both conservative and simple to 
interpret. Given real data, it would be worth putting substantial effort into understanding 
how detector and physics effects affect the shape of each distribution in turn. This could 
significantly improve particle mass resolutions. The resolutions obtained from the fits to 
the distributions in Figure ^ and Figure ^ are listed in Table ^. 

The Mt2 distributions obtained after the cuts listed in Section 3.4.1 are shown in 
Figure ^. The signal dislepton region having the desired edge is the unhatched region on 
each plot. Conveniently, most of the other SUSY events passing the cuts (mainly events 
containing gauginos) are also distributed with an edge located at a similar position to the 
disleptons'. This is due to the fact that such events often include pairs of slepton decays, 
and these may sometimes occur in combination with low jet activity and without additional 
leptons being produced inside detector acceptance. Consequently, for much of the SUSY 
background which passes cuts, the cuts are accepting Mt2 values which do not degrade 
edge performance significantly. Note that AM(x) = MT2ix) — X ^ is plotted evaluated 
at X = ' ^° should have an edge at AMmax = f^i — . Of course, m^o is not 
actually known a priori, so generating this graph in a real experiment requires an estimate 
of m^o to be obtained by other means. For our purposes it will be sufficient if the position 
of the edge of the distribution can be measured to about 10%, and this determines the 
accuracy required of the neutralino mass estimate x- The width AM{x) = Mt2{x) ~ X 
remains approximately stable at the 10% level for < x < 2m^o. To illustrate the 

relative insensitivity of AMmax(x) to X near rn^o, plots generated from the same data as 
before but with x ~ it^^o ± 50 GeV are shown in Figure ^. Satisfying the 10% requirement 
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above by obtaining a suitable value of x within such a ±40% range is not difficult and 
may be accomplished by first performing a "cut down" version of the analysis described 
later in the text, but omitting the Mt2 data, and then choosing the estimate, X; to be the 
resulting reconstructed neutralino mass. 

Determining the likely resolution for the Mt2 edge requires a slightly different ap- 
proach to that used for the other edges. Whereas the cuts used to obtain the invariant 
mass distributions have high SM rejections (primarily due to the presence of at least one 
high mass "pT or Meffective cut), the Mt2 cuts cannot be so hard, primarily because the 
desired dislepton events have have very little jet activity. With the dislepton production 
cross sections typically two orders of magnitude smaller than the squark/gluino production 
cross sections (Table |2|) a low efficiency is not affordable. There are also irreducible SM 
backgrounds (primarily T4^"'"Ty~ — > but also — > bbW~^W~ jjl'^l~vD in cases 

where jets are below the pT cut or outside detector acceptance) which have signatures 
identical to dislepton events. These backgrounds are clearly visible in Figure |^ and would 
cause problems for naive straight line fitting techniques. 

The standard model backgrounds in Fig- 
ure ^, although large, are easy to control. The 
SM edge is very clean, since the events which 
pass the cuts are well reconstructed and there 
are no noticeable tails. The SM edge will fall 
at mw for X = 0) corresponding to the mass 
of the missing neutrino in SM events. As x iii" 
creases, the SM edge recedes to lower masses 
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as shown in Figure 10, falling to 60 GeV for 
X = 100 GeV. A significant excess of events 
above this threshold would be a clear signal for 
non-SM processes. Such an excess will appear 
when the /-x? mass difference exceeds 80 GeV 
for low m^o and 50 GeV for high m^o. 

We estimate that, using Mt2 with either 
set of cuts, it is possible to measure AMmax = 
mi - m^o to 10% or better. The "hard" Mt2 
cuts successfully remove almost all SUSY back- 
ground above the SM threshold at the expense 
of only retaining half of the events passing the 
"soft" cuts. For both sets of cuts the SM thresh- 
old, at about 60 GeV, would be approximately three sigma away from AM^ 
accuracy. 
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Figure 10: Variation with x of the SM con- 
tributions to max(MT2(x) ~ x) (for a number 
of events corresponding to 100 fb~^) which 
provides a measure of the minimum l-Xi mass 
difference which is needed for signal Mt2 
events to be able to extend beyond the SM 
backgrounds. Events contributing to this plot 
were selected using the "hard" Mt2 cuts. 
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Figure 8: These Mt2 distributions for 01 are generated from 100 events fb ^ at high luminosity, 



and from the Mt2 cuts described in section B.4. (a) and (b) use the "hard" cuts, while (c) and (d) 
use the "soft" cuts. In (a) and (c), events containing a jet with pT > 40 GeV were vetoed. In (b) 
and (d), this cut was relaxed to 50 GeV. In all plots, Mt2{x) is evaluated at x = Ti-no. The results 



are presented in the form MT2{'mx°) 



in order to show the edge of the signal region located 



at the difference between the slepton and ncutralino masses (92.5 GeV). 
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GeV) is not presupposed. Instead, (a) uses x = 70 GeV and (b) uses x = 170 GeV. 



Endpoint 



edge 
l~^l~q edge 
l^q high-edge 
l^q low-edge 
l^l^q threshold 

Xq edge 
AM {Mt2 edge) 



S5 



Fit 


Fit error 


Fit 


109.10 


0.13 


70.47 


532.1 


3.2 


544.1 


483.5 


1.8 


515.8 


321.5 


2.3 


249.8 


266.0 


6.4 


182.2 


514.1 


6.6 


525.5 



01 

Fit error 



0.15 
4.0 
7.0 
1.5 
13.5 
4.8 
10% 



Table ^ values 
S5 01 



109.12 
536.7 
464.2 
337.0 
264.9 
509.2 
35.7 



70.50 
530.5 
513.6 
231.3 
168.1 
503.4 
92.5 



Table 5: Endpoints and associated fit uncertainties. The results of the naive fits shown in Figures ^ 
and ^ may be seen in the 'fit' and (one-sigma) 'fit error' columns. For the reasons outlined in 
Section 3.1 there is presently a lack of good theoretical predictions for these edge positions. (In this 
context a "good prediction" is one capable of taking into account, up to a level compatible with 
the statistical/fit errors, either the naivety of the fits, or the effects introduced by edge distortions 
and the presence of backgrounds.) The best guides presently available are the quantities listed 
in Table evaluated at rriq equal to the largest squark mass. These values are also shown for 
comparison purposes. All values are GeV except where stated otherwise. 
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3.6 Reconstructing sparticle masses 



Sparticle masses are reconstructed by performing a chi-squared fit with between four and 
six free parameters p. The main parameters of the fit are th^q, m^, m^o and rriq. Where 
appropriate, mh and mz also appear as fit parameters, although if they do, they are 
strongly constrained (particularly in the case of the Z) by present LEP measurements or 
expected LHC errors (0.0077% and 2.0% respectively). 

The chi-squared, as a function of the free parameters p, is then formulated as: 



ior-o,{p)y 



(3.8) 



where Of™ = Oi(pinodei) + criXi are n "smeared observables" , Xi ~ A^(0, 1) are the n 
random variables from the Normal distribution (with mean and standard deviation 1) 
which accomplish the smearing, ai is the anticipated statistical error for the i*^ observable 
(approximated by the fit error listed in Table ^), and Oi{p) is the value one would expect 
for the i**^ observable given the pcirameters p (for examples see Table Pmodei 

indicates 

the actual masses of the sparticles in the particular model being considered. The results 
of the fit, pfit, are then those which minimise the chi-squared. pfn may be interpreted as 
the sparticle masses which might be reconstructed by an LHC experiment after obtaining 
lOOfb^^ at high luminosity, with the Xi parametrising the experimental errors. In order to 
determine the accuracy with which this reconstruction can be performed, the above fitting 
process is repeated many times for different values of the Xi, producing the distributions 
of reconstructed particle masses which follow. 
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0.034 
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2.4 


22 GeV 
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Table 6: RMS and mean values, for each of the particles p, of the fractional mass-error distributions 



{Sp) in Figures 13 and |l4 Also shown are the widths of the distributions of the reconstructed particle 



masses (from Figures |ll| and . 

Figures ^ and |l^ show the probability distributions expected for the reconstructed Xi; 
/, X2 q masses, while Figures |l^ and |lj show the corresponding fractional errors, Sp, 
for the same quantities. All are approximately Gaussian and have reassuringly small tails. 
Statistics summarising these plots (widths, means and RMS values) are listed in Table ^. 

It may be seen that in both 01 and S5 the widths of the mass distributions for all 
four particles are very similar. This is because, in localized regions of parameter space, 
the edges tend to constrain mass differences far better than absolute masses. Evidence of 



this may be seen in Figure 15 which shows the scatter of reconstructions in the m^o-mj^ 
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Figure 11: Reconstructed Xij h X2 ^^^d 9 masses at S5. The small arrows indicate the masses 
used as the input parameters, pmodoi- 

plane. It is interesting to note that without the Mt2 constraint at 01, the fit's chi-squared 
commonly has two distinct and competing comparable minima - one at high and one at 
low values of tni^. The simultaneous existence of these minima is a direct consequence 
of putting the near and far l^q edges on an equal footing in this analysis, allowing more 
than one interpretation for each of the high and low edges. The need to resolve this kind 
of ambiguity in model-independent investigations of this type illustrates the importance of 
establishing reliable model-independent ways of measuring the absolute scale of m^^ — m^o 
(or a related quantity) even if only to an accuracy of 20-30%. 

The clear gap between the reconstruction regions for 01 and S5 in Figure ^ supports 
the original claim that, systematic errors permitting, it will be possible to distinguish 
between the two scenarios. 
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Figure 12: Reconstructed Xii X2 ^-^^d q masses at 01. The small arrows indicate the masses 
used as the input parameters, Pmodci- 
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Figure 15: Reconstructed mj versus m^o for 01 and S5. The stars show the true masses for each 
modeL 
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4. Conclusions 

Some of the five standard LHC SUGRA points are compatible witli universal perturbative 
string and M-theory, but dangerous CCB/UFB breaking minima are present in each exam- 
ple. We therefore studied a perturbative string model which is optimized to ameliorate the 
CCB/UFB problems present in the other models. The optimized model is non-universal 
because the squarks and sleptons are split in mass at the string scale. We identify the 
SUGRA point with the most similar spectrum and hard SUSY production cross sections 
(S5) to compare the optimized model with. The main difference is that the sleptons are 
heavier and therefore have lower production cross-sections. We have demonstrated the 
existence of a method by which an LHC experiment will be able to measure the masses of 
the (lighter) sleptons and the two neutralinos at 01 in a largely model independent way. 
In a specific comparison of S5 and 01 we have shown that this method will be able to 
distinguish a SUGRA model from an optimised string model with very similar properties. 

More importantly, we expect that the techniques developed here are general enough 
to be used to discriminate between other pairings of optimised and non-optimised models 
with similar characteristics. The optimized model analysis applies to a more general class 
of models than the string model itself. We could apply it to models with non-universal 
SUSY breaking terms at Mqut in which the squarks and sleptons are explicitly split in 
mass. These constitute a superset of the particular string model considered here. 
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